Peeling Bifurcations of Toroidal Chaotic Attractors 
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Chaotic attractors with toroidal topology (van der Pol attractor) have counterparts with symme- 
try that exhibit unfamiliar phenomena. We investigate double covers of toroidal attractors, discuss 
changes in their morphology under correlated peeling bifurcations, describe their topological struc- 
tures and the changes undergone as a symmetry axis crosses the original attractor, and indicate 
how the symbol name of a trajectory in the original lifts to one in the cover. Covering orbits are 
described using a powerful synthesis of kneading theory with refinements of the circle map. These 
methods are applied to a simple version of the van der Pol oscillator. 



I. INTRODUCTION 

It has been known for some time that discrete symme- 
try groups can be used to relate chaotic attractors with 
different global topological structures [l|, [2, H, 0] . By dif- 
ferent (or distinct) topological structures we mean there 
is no smooth deformation of the phase space that can 
be used to transform one attractor into the other in a 
smooth way. If a chaotic attractor has a discrete sym- 
metry, points in the attractor that are mapped into each 
other under the discrete symmetry can be identified with 
a single point in an "image" attractor. The identifica- 
tions are through local diffeomorphisms. The original 
symmetric attractor and its image are not globally topo- 
logically equivalent. This process can be run in reverse. 
A chaotic attractor without symmetry can be "lifted" to 
a covering attractor with a discrete symmetry following 
algorithmic procedures [H, 0, H, 0, [1] . 

A simple example illustrates these ideas. The Lorenz 
attractor [6| obtained with standard control parameter 
values exhibits a two-fold symmetry. The symmetry is 
generated by rotations about the Z axis through tt ra- 
dians: i?z(7r). We mod out this two-fold symmetry by 
identifying pairs of points {X, Y, Z) and {—X, —Y, Z) in 
the symmetric attractor with a single point (u, v, w) = 
{X^ — Y^, 2XY, Z) in the image attractor. This results in 
a chaotic attractor that is not topologically equivalent to 
the original attractor. Rather, it is topologically equiva- 
lent (not diffeomorphic, [H, 0, 0]) with the Rossler attrac- 
tor [3]. Similarly, the Rossler attractor can be lifted to 
a two-fold covering attractor that is topologically equiv- 
alent to the Lorenz attractor. The lift is carried out by 
inverting the 2-^1 local diffeomorphism used to gen- 
erate the image: {u,v,w) {X = ±^^{r + u), Y = 



Z 



where 



X^ +Y^ 



This modding out process is illustrated in Fig. [TJ 

A single image attractor can have many topologically 
inequivalent covers, all with the same symmetrygroup. 
These covers are differentiated by an index 0, [J- The 
index has interpretations at the topological, algebraic, 
and group theoretical levels. Briefly, the index describes 



FIG. 1: The Lorenz attractor (a) can be mapped to a 
Rossler-like attractor (b) by identifying points related by ro- 
tation symmetry about the Z-axis. This process is reversible: 
Rossler-like attractors can be "lifted" to Lorenz-like attrac- 
tors. 



how the singular set of the local diffeomorphism relating 
cover and image attractors is situated with respect to the 
image attractor. Different lifts of the Rossler attractor, 
all with Rzi'K) symmetry but with different indices, are 
obtained if the two-fold rotation axis: passes through 
the hole in the middle of the Rossler attractor; passes 
through the attractor itself; or passes outside both the 
attractor and the hole in the middle (cf., Fig. The 
transition of the symmetry axis through the attractor is 
responsible for peeling bifurcations 0]- 

Rossler-like attractors have been lifted to covers with 
many symmetry groups Whenever the singular set 
of the symmetry group involves a rotation axis, this axis 
has passed through the attractor at most once in all pre- 
vious studies. More precisel y, it has passed through the 
branched manifold [1, [lol. Illl| describing the attractor 
at a single point. 

There have been no studies of covering attractors that 
are obtained when the rotation axis intersects the image 
attractor in more than one spot. Generically, a rotation 
axis must intersect a toroidal attractor an even number 
of times. In the present work we look at two-fold covers 
of toroidal attractors with rotation symmetry. We use 
methods similar to those used in 111 vA, i S H. Many of 
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FIG. 2: The Rossler attractor (a) can be lifted to topologically 
distinct double covers with rotation symmetry by placing the 
rotation axis in different positions: (b) in the attractor; (c) in 
the hole in the middle of the attractor. 



our results depend on a powerful synthesis of kneading 
theory with refinements of the circle map. 

In Sect. IVII below we study peeling bifurcations for 
rotation-symmetric lifts of a chaotic attractor generated 
by the van der Pol equations subject to a periodic drive. 
The phase space of this attractor is a "hollow donut" or 
"fat tire" , that is, the direct product of an annulus with 
a circle: x S^. An annulus itself is a circle with a 
thick circumference; = x /, where / is an interval, 
or short segment of the real line. The branched manifold 
that describes this attractor is essentially a torus (thin 
tire). A rotation axis used to construct symmetric lifts 
must intersect the torus an even number of times. Lifts 
of laminar and chaotic flows on a torus are discussed in 
Sees. HV] and El respectively. The first of these sections 
describes how the symbolic dynamics of image trajec- 
tories lift to symbolic dynamics in covering trajectories. 
The second describes how the branched manifold describ- 
ing the image attractor lifts to the branched manifold 
describing the covering attractor. We prepare for these 
discussions by introducing toroidal coordinates and de- 
scribing flows on a torus in Sec. IIIII We begin this entire 
odyssey in Sec. with a review of peeling bifurcations 
and their properties. Our results are summarized in Sec. 

Em 



II. REVIEW OF PEELING BIFURCATIONS 

Peeling bifurcations arise naturally when considering 
covers of chaotic attractors. They describe the bifurca- 
tions these covers can undergo as the relative position 
of the image attractor and the symmetry axis changes. 
Peeling bifurcations have been described in some detail 
for covers of the Rossler dynamical system in [1, [J ■ 

We briefly describe the basic idea for double covers 
with i?z(7r) symmetry about a rotation axis TZ with the 
usual saddle-type symmetry. Lift an image attractor 
(Fig. [2Ia)) to a double cover with TZ far away from the 
original image attractor. The double cover consists of two 
identical copies of the original image attractor. They are 
disconnected. An initial condition in one will evolve on 
that attractor for all future times in the absence of noise. 
As the rotation axis TZ approaches the image the two 
disjoint components of the double cover approach each 
other, keeping TZ between them. At some point TZ will 
intersect the attractor. When this occurs the rotation 
axis will split the outer edge of the flow from one of the 
two components of the cover and send it to the other com- 
ponent, and vice versa. The two attractors in the cover 
are no longer disconnected (Fig. [2Ib)). As the rotation 
axis TZ moves deeper to the center of the image attractor 
(towards the center of rotation of the Rossler attractor, 
for example), the double cover becomes smaller in spatial 
extent. Finally, the rotation axis TZ may stop intersect- 
ing the image by passing into the hole in the middle (Fig. 
He)). 

The peeling bifurcation takes place as the rotation axis 
moves from the outside to the inside of the image attrac- 
tor. The image attractor itself is not affected. All bifur- 
cations take place in the cover. Before intersections be- 
gin and after they end the double covers are structurally 
stable and topologically inequivalent. During the inter- 
section phase the cover is structurally unstable because 
slight changes in the position of the rotation axis produce 
profound changes in lifts of trajectories from the image. 
That is, a trajectory remains unchanged in the image 
while its lift(s) into the cover change dramatically as the 
rotation axis TZ moves. Changes in the shape, structure, 
and period of lifts of unstable periodic orbits from the 
image into the cover are predictable. These changes are 
summarized for covers of Rossler-like attractors in Figs. 
7 - 12 of ref. 0. 

It may be useful to regard peeling bifurcations in terms 
of how two trajectories in the cover that result from the 
lift of a single trajectory in the image connect or re- 
connect as the position of the symmetry axis changes. 
Roughly but accurately, they can 'turn back' into the 
subset of the attractor from which they originated when 
the axis is outside the trajectory, or else 'cross over' into 
the complement of that subset when the axis is inside the 
trajectory. This behavior is reminiscent of what happens 
during the transition of a plane through a saddle point 
on a surface, with the added feature of 'direction'. Since 
orbits are dense in a strange attractor, the lifted system 
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cannot be structurally stable during a peeling bifurca- 
tion. 

III. FLOWS ON A TORUS 



done by setting p — 2,r = \,6 = acj) in the toroidal coor- 
dinates. This curve closes or does not close depending on 
whether a is rational or irrational. The return map on a 
plane cj) = const, is shown in Fig. [H It is 6n+i = 9n + n 
mod 27r. 



It is useful to describe flows on a torus in terms of a 
system of coordinates adapted to the torus: {(j),9,r). In 
such a coordinate system (j) is the longitude; it increases 
with time: d(j)/dt > 0. The angle 9 is the meridional 
angle, measured from "the inside of the torus" (see Fig. 
[3]), and r measures the distance from the center line of 
the torus, a circle of radius p in the x-y plane. The circle 
radius p must be sufficiently large so that p — r>0 for all 
points in the attractor. Standard cartesian coordinates 
are represented in terms of these toroidal coordinates by 

X — {p — r cos 9) cos 

y = {p — r cos 9) sill (f) (1) 

z = — rsin^ 




/2k 

n 

FIG. 4: The return map for a rigid rotation is a straight line 
mod 2tt. Q./2-K = 0.20. 




FIG. 3: Toroidal coordinates. 



The Birman- Williams theorem [1, [T^, [Tl[ can be 
applied to dissipative toroidal flows in that generate 
strange attractors. The result is that the topology of the 
flow is described by a branched manifold. The mecha- 
nism generating chaotic behavior involves an even num- 
ber of folds. The branch "lines" are now circles. Since 
the flow occurs in a bounding torus [13,1131 genus one, 
the Poincare surface of section consists of a single disk. 
The intersection of the branched manifold with the disk 
("branch fine") is topologically a circle, . As a result, 
the flow can be investigated by studying maps of the cir- 
cle to itself [Hi. 

IV. LIFTS OF RIGID ROTATIONS 

In order to determine the topological structure of a 
strange attractor in it is sufficient to determine the 
topological structure of the branched manifold that de- 
scribes it. This remains true for covers of strange attrac- 
tors with arbitrary symmetry H, fl, 13] . We do this in the 
following section. 

In this section we prepare the groundwork by inves- 
tigating how a rigid rotational (quasipcriodic) flow on a 
torus is lifted to a double cover of the torus. This is easily 



Now pass a rotation axis through the torus as shown in 
Fig. O The order-two rotation axis intersects the torus 
in an interval /. If the rotation axis is parallel to the Z 
axis, the end points of this interval are at 6* = 27r(i ± a). 
The rotation symmetry lifts the torus into a structure 
inside a genus-three bounding torus that is shown in Fig. 
[6] The location of the rotation axis is indicated by x . 




27i(l/2-o) 



FIG. 5: The order- two rotation axis intersects the torus at 
27r(i ±cr). 




FIG. 6: When the order-two rotation axis intersects the torus 
at 27r(i ± a), the double cover consists of a structure with 
the topology shown. This geometric structure exists inside a 
bounding torus of genus three. 

Since the flow exists in a bounding torus of genus three 
the global Poincare surface of section has two components 
Iisj . In such cases the first return map consists of a 
2x2 array of maps [1, [1^1 ■ For the Lorenz attractor such 
maps indicate flows from branch line to branch line. In 
the present case the return map indicates flows from the 
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FIG. 7: The laminar flow on a torus lifts to a flow on a man- 
ifold with the complicated form shown. The rotation axis is 
indicated by cross. 



branch circles on the left and right of Fig. [71 The return 
map for the cover of the rigid rotational flow, in the case 
that the Z axis intersects the image torus at 9 = 7r±27rt7, 
is presented in Fig. [H The angles parameterizing the 
branch circles on the left and right run from zero to 27r. 
This figure shows that an initial condition within 27rcr of 
9 = n on the left hand branch circle maps to the right 
hand branch circle (see Fig. [5]Ja)), and vice versa. 

A symbol name for any trajectory on the covering flow 
is easily constructed. Assume the return map in the im- 
age is 9n+i = 9n + ^ mod 2tt. Then 9n ^ 9a + nil mod 
2TT. Write out this string of real numbers and replace 
each value 9n by 1 ii 9n & I , zero otherwise. This results 
in a string of symbols, for example 00000 11111 00000... 
for initial condition L. Choose an initial condition L or 
R, for one side of the cover or the other. Then repeat 
this letter following symbol 0, conjugate this letter (L 
i?, i? — *■ L) following symbol 1. This algorithm leads 
to 000001111100000.. LLLLLRLRLRRRRRR... A 
rotation-symmetric trajectory has a conjugate sequence. 

Depending on parameter values (e.g., a <C tr < j) the 
symbol sequence can consist of long strings of is, long 
strings of i?s, and long strings of LRs, giving the appear- 
ance of prolonged rotation about three centers: one being 
the left-hand torus in the lift, another being the right- 
hand torus in the lift, and the third alternation about 
both in sequence when 0„ falls in the interval / over a 
large range of successive interations. 




L R 



FIG. 8: Return map for a rigid flow contained within a genus- 
three torus consists of mappings from two branch circles to 
themselves. 



As the rotation axis sweeps from the outside to the 
inside of the torus, the value of a increases. The circu- 
lar intervals for which transition from one side to the 
other takes place increases, with the return map be- 
coming more and more off-diagonal. As the Z axis ap- 
proaches the inner part of the image torus, the measure 
of 9 values that map to the same branch circle decreases, 
and becomes zero when tangency occurs (ct = tt). At this 
point the return map is completely off-diagonal. This is 
an indication that the global Poincare surface of section 
is no longer the union of two disjoint disks. A single disk 
suffices. This signals that the flow returns to a flow of 
genus-one type, and the return map on the single disk is 
9n+i = 9n + 2Q mod 2tt (notice the factor of 2). 

In the limit when the rotation axis is outside the torus 
("cr < 0") the double cover consists of two disconnected 
tori. An initial condition in one torus remains forever 
in that torus. In the limit when the rotation axis is in 
the hole in the middle of the torus ("cr > i") the double 
cover consists of a single torus. When the rotation axis 
goes through the origin of cartesian coordinates, the lon- 
gitudinal angle (f> in the image increases twice as fast as 
the longitudinal angle $ in the cover. 

Simulations of peeling bifurcations for double- and 
triple-covers of laminar flows on a torus can be found 
at 0. 



V. LIFTS OF CHAOTIC ATTRACTORS 

There is a class of strange attractors, such as the van 
der Pol attractor that we discuss in the following sec- 
tion, whose phase space is a hollow donut, topologically 
X S^, where describes the longitudinal (flow) di- 
rection. The intersection of the attractor with a constant 
phase plane (j> —est. occurs in an annulus = x I, 
where this describes the meridional direction and / is 
a small interval. Under the Birman- Williams projection 
fldl . [Tl| the intersection of the projected attractor with 
the plane 4> =cst. is topologically a circle S^. The for- 
ward time map (j> <j) + 2tt is therefore a mapping of 
the circle to itself. As a result, many of the properties of 
this class of attractors are determined by the properties 
of the circle map. For example, the number and labeling 
of the branches of the attractor's branched manifold are 
determined by the appropriate circle map. 

For simplicity we assume that three branches A, C, B 
suffice to describe the branched manifold, and that the 
return map on the singularity at which these branches 
are joined (the branch circle) is a circle map [l3 |: 

9n+i = 9n + ^ + Ksm9n mod 27r (2) 

with K > I. This map is shown in Fig. [51 Branch 
C is orientation reversing. Its forward image extends 
over a range less than 2tt and its extent is delineated by 
the two critical points. Branches A and B are orienta- 
tion preserving and their forward image extends over a 
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range greater than 27r. They are dehneated by the criti- 
cal points that bound C and the inflection point between 
them. 
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FIG. 9: Circle map Eq. ^ ioT K = 2.5, n/27r 
three branches are conveniently labeled A,C,B. 



0.2. The 



The return map for the double cover is obtained as in 
the previous section. We begin by looking at intersections 
of the rotation axis near the outside of the torus, at values 
9 = TT ± 27rcr, with a small. The return map on the two 
branch circles is as shown in Fig. [10] for a — 0.15. 













/ 
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FIG. 10: Return map for the double cover of the chaotic 
flow whose return map is shown in Fig. (5] Parameter value: 
a = 0.15. 



The return map for the double cover is obtained from 
the return map for the image as follows. The vertical 
lines through the maximum and the minimim and the 
vertical axis at 0/27r = 0, 1 in Fig. [9] separate the return 
map into three branches A, B, C. 

The two additional vertical lines separate branch C 
into three branches: C3 which is the interval /: 7r — 27r cr < 
9 < IT + 27rcr, and Ci and C3, which map L ^ L and 
R ^ R. Generally there are no degeneracies, so these 
five vertical lines divide the circle into five angular inter- 
vals. In the present case both endpoints of the interval / 
occur inside the orientation-reversing branch C, so this 
branch is divided into three parts: Ci, C2, C3, as shown 
in Fig. [9l As a result, initial conditions on branches A 
and -B, and the adjacent parts of branch C, namely Ci 
and C3 of the circle on the left, map back to that circle 



while initial conditions in the angular interval C2 on the 
left circle map to the right hand circle. The branched 
manifold describing the covering flow has 10 branches 
with transitions summarized as follows: 

L R^ R L-> R R~^ L 



A 

C3 
B 



A 

C3 
B 



Co 



In the event that one (both) of the endpoints of the in- 
terval / coincide with one (two) of the three points sep- 
arating A, C, B there are 8 (6) branches. 

Before the peeling bifurcation begins, when the two 
identical covers are well separated, each is charac- 
terized by a branched manifold with three branches. 
At the end of the peeling bifurcation, when the ro- 
tation axis 'R- is inside the image torus, there are 
9 = 3^ branches that can be labeled (^4, C, B) (g) 
{A,C,B) = [AA,BA,CA,AB,BB,CB,CA,CB,CC] 
[3, Branches labeled by an even number of letters 
C are orientation preserving. 

A periodic orbit in the image can be lifted to one or 
two covering orbits. The symbol name of the covering 
orbit is obtained from the symbol name of the image or- 
bit and information about the interval /. The name of 
the orbit in the image is written out (e.g., ABBCBBAC 
and refined according to the location of the interval / 
(e.g., ABBC2BBAC3. An initial condition on the left or 
right {L or R) is given and this is changed whenever a 
trajectory passes through one of the branches defined by 
/. For example, under this algorithm ABBCBBAC 
AlBlBlC2lBlBlAlC3l ARB]iBjiC2RBRBjiAjiC3R. 
A period-p orbit lifts to two period-p orbits or one sym- 
metric orbit of period 2p depending on whether the image 
orbit maps through the interval / an even or odd number 
of times. The algorithm for lifting orbits from a toroidal 
fiow to a double cover involves a synthesis of kneading 
theory with refinement of the circle map due to the in- 
tersection of the rotation axis with the image toroidal 
flow. 



VI. APPLICATION TO THE VAN DER POL 
ATTRACTOR 

The Shaw version version of the periodically driven 
van der Pol equations 



u — bv + {c— dv'^)x 
V — —u + asm{ujt) 



(3) 



produce a toroidal attractor for control parameter values 
(a, 6, c, d, uj) = (0.25, 0.7, 1.0, 10.0, 7r/2) [1, Q. The phase 
space for this attractor is the direct product of an annular 
disk with a circle. One projection of this attractor is 
shown in Fig. [TlJ 
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u 

FIG. 11: Projections of the van der Pol attractor 
on the plane u-v. Parameter values: (a, 6, c, d, u) = 
(0.25,0.7, 1.0, 10.0, 7r/2). 



The attractor is mapped into R'^ following the pre- 
scription x{t) — {p — u{t)) cos{ujt),y{t) = {p ~ 
u{t)) sm{ujt), z{t) = v{t), with p — 1.2. This flow, em- 
bedded in K^, has the topology of a hollow donut. A 
projection onto the x-y plane is shown in Fig. 1121 




-2 -1 1 2 

X 

FIG. 12: Mapping of the van der Pol attractor into using 
Eq. dl} with p = 1.2. 

The covers of the chaotic attractor produced by this 
embedding into undergo correlated peeling bifurca- 
tions as the rotation axis slices through the image. In 
Fig. [T3]we show three projections of the double cover ob- 
tained when the symmetry axis is parallel to the Z axis 
and has two-fold symmetry. The rotation axis passes 
through the point {x,y) — (1.20,0.0) in the x-y plane. 
The continuous version of this correlated peeling bifur- 
cation is available at |16l] . 

Both this attractor and the Lorenz attractor are con- 
tained in genus-3 bounding tori. The two attractors are 
topologically incquivalcnt. The X-Y projection of this 
attractor, which is shown in Fig. [TSja), is similar to 
the X-Y projection of the Lorenz attractor. However, 
projections onto the other two directions are totally dif- 
ferent. The toroidal structure of the present attractor is 
revealed in the projections shown in Figs. [TSKb) and (c). 
A Poincare section of the double cover of the van der Pol 
attractor (Fig. [H)) shows the double annular shape. The 
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FIG. 13: Double cover of the chaotic attractor solution of 
Equ.Q for the Shaw version of the van der Pol equations is 
mapped from xS^ by a natural embedding. The center line 
of the torus is mapped to a circle of radius 1.2 in the x-y plane. 
A correlated peeling bifurcation occurs when the double cover 
is around the two-fold rotation axis through (x, y) = (1.2, 0.0). 
The three images are projections onto the X-Y, X-Z, and Y- 
Z planes, from top to bottom. 
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two components of the Poincare surface of section consist 
of two half planes, both with y = 0. One is hinged on 
an axis parallel to the Z axis through [X^Y) = (1.1,0); 
the other is the rotation image of the first. Intersections 
with y = 0, y > are taken on one half-plane and in- 
tersections with y = 0, y < are taken with the other. 
This Poincare section emphasizes the invariance of this 
attractor under rotation symmetry around the Z axis. 



and those of its covers. Specifically, an attractor whose 
phase space is a hollow donut intersects a rotation axis 
an even number of times (more precisely, its branched 
manifold does). Further, its branched manifold can have 
only an odd number of branches. These remarkable prop- 
erties extend, in a suitable way, to double covers of these 
attractors. 




FIG. 14: Intersections of the double cover (Fig. I13|l of the 
van der Pol attractor with the two disconnected components 
of the Poincare section. The component on the right has 
y = 0,X > +1.1 with Y > and the component on the left 
has y = 0, X < —1.1 with Y < 0. The rotation symmetry is 
clear. 



VII. SUMMARY 



For the first time methods for constructing double cov- 
ers of chaotic attractors have been applied to chaotic at- 
tractors of a toroidal nature. These attractors are con- 
tained in genus-one bounding tori and are described by 
branched manifolds with a circular cross section on a 
Poincare surface of section [9|. Their return maps are 
maps of the circle to itself. Their double covers are cre- 
ated by correlated peeling bifurcations. The morphology 
of the covering attractor changes systematically as the ro- 
tation symmetry axis slices through the image torus from 
outside to inside. Outside, the double cover consists of 
two identical attractors, each contained in a genus-one 
torus. The two genus-one tori are disconnected. When 
the rotation axis intersects the image, the double cover is 
contained in a genus-three torus and is not structurally 
stable against perturbations of the position of the ro- 
tation axis. When the rotation axis enters the hole in 
the torus, the double cover exists in a genus-one torus. 
For various ranges of lift parameter values rotations can 
appear to occur around a single center, two centers, or 
three centers. Lifts of periodic orbits in the image at- 
tractor are described by a powerful synthesis of kneading 
theory with refinements of the circle map. 



It is remarkable that the global topology of the image 
attractor imposes nontrivial constraints on its properties 
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